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Abstract 

There is an observational indication of extragalactic magnetic fields. No known 
astrophysical process can explain the origin of such large scale magnetic fields, which 
motivates us to look for their origin in primordial inflation. By solving the linearized 
Einstein equations, we study metric perturbations sourced by magnetic fields that are 
produced during inflation. This leads to a simple but robust bound on the inflation 
models by requiring that the induced metric perturbation should not exceed the 
observed value 10 -5 . In case of the standard single field inflation model, the bound 
can be converted into a lower bound on the Hubble parameter during inflation. 



1 Introduction 



It has been well known that magnetic fields of ~ /iG exist in galaxies (see for example, 
P-H]). These fields are thought to be an outcome of the amplification of seed magnetic 
fields of unknown nature through the dynamo mechanism. Magnetic fields are also known 
to exist in cluster of galaxies, too pQ. Although such seed magnetic fields are expected 
to still reside in the intergalactic medium and in voids, there had been no observational 
support of the seed magnetic fields until quite recently. In 2010, observational data taken 
by Fermi and High Energy Stereoscopic System (HESS) gamma-ray telescopes provided 
a strong support of the existence of extragalactic magnetic fields of at least B ~ 10~ 17 G 
on Mpc scales [5H7]. The lower bound on B was derived from detection by HESS of TeV 
gamma rays coming from TeV blazars and non-detection by Fermi of GeV scale cascade 
emission, which is in contradiction with the assumption of zero magnetic field. 

Explaining the origin of the extragalactic magnetic fields of B ~ 10~ 17 G on Mpc 
scales is challenging and remains to be addressed. No promising astrophysical processes 
are known to generate the suggested amplitude of magnetic fields on such large length 



scales Ej. In light of this situation, the best alternative we can think of is to make use of 
inflation in the early Universe [T8H20] which causally connects length scales far beyond the 
Hubble radius through superluminal expansion [2"TH2"9"] . It is well known that to realize 
magnetogenesis by inflation the action for the electromagnetic field needs to be modified 
such that the action breaks conformal invariance. The most-used Lagrangian for this 
purpose is a form / 2 (0)F /ii/ F /ii/ , where <fi is some scalar field (e. g. an inflaton or a 
dilaton) that varies during inflation. It is not easy, however, to obtain a successful model, 
because typical models encounter either breakdown of perturbation theory due to the 
strong coupling that occurs in the earlier stage of inflation [30], or generation of excessive 
electric fields whose energy density may cause serious backreaction [2T||23]. 

The purpose of this paper is to discuss the latter effect caused by primordial magnetic 
fields without respect to its specific generation mechanisms. As already discussed in the 
literatures [2"9" |l3"TH3"o*] . if the energy density of the magnetic field exceeds the background 
inflaton energy density, its backreaction onto the background dynamics becomes signifi- 
cant, which generally destroys the homogeneity and isotropy of the Universe and puts a 
constraint on model building. In addition to this, even if the backreaction effect is small 
and is not problematic for background evolution, there is another issue that one must take 
into account. The energy-momentum tensor of the magnetic field also induces metric per- 
turbations and distorts the homogeneous and isotropic FLRW Universe [3"oTl3"%] . Degree 
of the distortion depends on the inflation models and also on the electromagnetic models 
that achieve magnetogenesis. But whatever be the model of inflationary magnetogene- 
sis, the resultant amplitude of metric perturbation must not exceed the observed value 
~ 10~ 5 [39]. Put it in another way, in addition to the standard constraints of explain- 
ing the observed feature of primordial perturbations, we can derive an additional generic 



^ tl For magnetic field generated from standard cosmological perturbations, see [BHTTj. For the effect of the 
primordial magnetic field on CMB temperature anisotropy and its observational constraints, see [T^HTT] . 
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bound that must be imposed on any inflation model by analyzing metric perturbations 
arising from the inflationary magnetic fields, which we work out in this paper. 

To achieve this in quantitative manner, we study evolution of metric perturbations 
induced by magnetic field which is supposed to be generated during inflation. We will 
treat the energy-momentum tensor of the magnetic field as first order perturbation and 
solve the linearized Einstein equations on the FLRW Universe throughout the inflation- 
ary era and the subsequent radiation dominated era. This paper reports the results of 
such a calculation and derives a new constraint on inflation models. It is found that the 
slow-roll parameter that characterizes how slowly the Hubble parameter changes during 
inflation needs to be bounded from below in order to avoid large metric perturbation from 
magnetic field. This is basically equivalent to the constraint 5p- m t > Sps which tells that 
the fluctuation of inflaton energy density must be larger than the energy density of the 
magnetic field. 

2 Metric perturbation from magnetic field 

In this section, we consider the evolution of the metric perturbation during inflation and 
subsequent radiation dominated Universe in the presence of magnetic field generated quan- 
tum mechanically by breaking the conformal invariance. Since magnetic field is assumed 
to be absent at the background level, we will treat its energy-momentum tensor as lin- 
ear perturbation. Perturbations of inflaton field and radiation/matter energy-momentum 
tensor are also treated as linear perturbations. 

Let us first write down the perturbed metric having only scalar-type perturbations in 
general gauge (for example, see [30]); 



ds 2 = -(1 + 2A)dt 2 + 2a 2 d i Bdx i dt + a 2 (t) [(1 + 2^)% + 2didjE] dx l dx j . (1) 



Here A, B, ip and E are metric perturbations. We denote the perturbed energy-momentum 
tensor including both matter such as inflaton or radiation and magnetic field as 



ST° = -Sp, 5T° l = {p + P)a 2 d l v, 5T\ = -{p + P)di(v - B), ST) = S^SP + n), (2) 



where v is defined by itj = a 2 d,iv and n l i = is the anisotropic stress. Since in many 
inflation models the anisotropic stress of matter is second order in perturbation, we will 
assume that only the magnetic field yields the first order anisotropic stress. If the time 
evolution of the magnetic field is solely due to the cosmological expansion, n % j decays 
in proportional to a~ 4 just in the same way as radiation. On the other hand, if the 
electromagnetic action is modified, for example, by having the electromagnetic field tensor 
couple to the inflaton like 




(3) 
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where <fi is inflaton, then can be made to decay much slower than a~ A . In particular, 
in order to have scale invariant power spectrum of magnetic field on cosmological scales, 
7T l j needs to be almost constant during inflation. This is indeed possible for some suitable 
choices of the functional form of /(</>). 

For later convenience, we will write down the explicit transformation properties of 
perturbation variables under the gauge transformation, t — > i = t — T, x % — > x l — x 1 — diL; 

A = A + f, B = B + L-^, i> = ij + HT, E = E + L, a g = a g -^, (4) 
. T 

5p = 5p + pT, 5P = 5P + PT, v = v-—, (5) 

or 

where a g = B — E denotes the shear of the four- velocity u M . 

The linearized Einstein equations in Fourier space are given by 

P 

—ip + Hk 2 a g - 3H 2 A + 3Hij; = 4nG5p, (6) 

(X 

HA-ij; = -4nG(p + P)a 2 v, (7) 
P P 



: (A + ip) + (3H 2 + 2H)A -i/j- 3HiP + HA (& g + 3Ha g ) = AnGdP, (8) 



3d 1 ' ' T ' ' ' 3 

P 

— (A + iP) + k 2 (a g + 3Ha g ) = SuGIl, (9) 



and n is defined by 



**, = ( 5«» " k ~W ) n - (10) 



The linearized conservation laws V M T^ = are given by 

S'p + 3H(6p + 5P) + 3(p + P)i; - (p + P)k 2 (v - <T g ) = 0, (11) 
d t [(p + P)a 2 v] +5P-\u + 3H{p + P)a 2 v + (p + P)A = 0. (12) 

To see how the curvature perturbation is affected by the presence of magnetic field, 
let us consider the curvature perturbation on the uniform energy density hyper-surface on 
which dp = is satisfied. Denoting the curvature perturbation on this slice by we can 
derive its evolution equation by using both Einstein equations and conservation laws given 
above. On superhorizon scales k/(aH) < 1, we find that the evolution equation reduces 
to 

1 / a 3 H Y 8ttG, 
a?\p + P C 



{ + 3H(+-[-— R 6P iel ) n = 0, (13) 



where SP rc \ = SP cm — jSp cm (Sp cm and SP cm is the energy density and pressure perturbation 
for the electromagnetic field, respectively) is the nonadiabatic pressure perturbation due 
to the relative entropy perturbation [4"0||4"T] between the magnetic field and the component 
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dominating the Universe [£3. Inhomogeneous solution of this differential equation is given 
by 



+ ~ 7I7 ^ a (*2)n(t 2 ) 



3 A. « 3 (^i) A, 



(14) 



where we have imposed an initial condition that C(^*) = 0- I n the actual situation, t* 
may be taken to be a horizon crossing time. Equation (I14p represents the magnetic field 
contribution to the curvature perturbation. 

Determining a precise value of the time integral in Eq. ( TT4|) requires both model specifi- 
cation and numerical computation, which is not what we want to do in this paper. Instead, 
without referring to the specific inflation model, we can estimate the order of magnitude 
on the basis of some generic properties of inflation. Let us first evaluate the first term 
containing SP Te \ in the radiation dominated Universe achieved after reheating. Noting that 
the electromagnetic tensor is traceless <5P re i during inflation can be approximated as 



5P r 



el — g^Pcm? 



(15) 



where we have used an approximate relation p ~ —P. In the radiation dominated Universe, 
the energy density of magnetic field decays in the same way as that of radiation and we 
have SP TC \ = 0. Therefore, the first term yields a constant: 



dti 



p(h) + P(h 



SPrelih 



2Af6 Pcm 



(16) 



where e = —H/H 2 is the slow-roll parameter, p in f is the energy density of the inflaton 
and M is the number of e-fold measured from the time when the mode of interest crossed 
the horizon to the end of inflation. Equation ( fl6l) becomes exact if all the quantities are 
constant during inflation but may deviate from the correct one if one of the quantities 
significantly changes during inflation. While it is a good approximation to treat p- m f and 
e as constants in many inflation models, this does not necessarily hold for 5p em . In order 
to be conservative as possible as we can, let us replace Eq. ([161) by 



H(h) 



dt-\ 

u p{t x )+P{t x 



SPrelih 



1 

> - 

e 



Sp em {t 



end J 



Pinf 



(17) 



# 2 Precisely speaking, intrinsic entropy perturbation for the dominant component Pintr will be correlated 
with dP r ei- But since the electromagnetic field is subdominant, it does not produce Pintr comparable to 
<5P rc i, or almost equal to — P re i- Therefore, we do not need to consider the intrinsic entropy perturbation 
for the dominant component. 

# 3 Traceless nature of the electromagnetic tensor persists even when the electromagnetic action is 
modified to a typical one given by Eq. ([3]). This may not be true for other modification of the Maxwell 
equations. Eq. (|15[) becomes an overestimate only if the effective equation of state parameter w m = P m / 1 p m 
for the electromagnetic field becomes very close to —1, which is unlikely. 
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where t en d is the time of inflation end. A merit of this replacement is that, as we will 
see later, |5p e m(^end)/Pinf| can be connected to the observed magnetic field strength and 
duration of the dust-like Universe after inflation but before the reheating. If Sp em is a 
growing function like 5p cm oc dP (p > 0), then left-hand side and right-hand side of the 
above inequality are almost the same. If it is a decreasing function, then left-hand side 
becomes much larger than the right-hand side. In either case, the left-hand side must 
not exceed the observed value of the primordial curvature perturbation which is about 
10 -5 . If we further use another observational constraint that the curvature perturbation 
is almost Gaussian [39]H2], the bound can be made by a few orders of magnitude tighter. 
This is because the curvature perturbation from Eq. ( TTTj) is proportional to the magnetic 
field squared and hence is highly non- Gaussian [29~ ] |4"3 ] |4*4"]. But since we do not know 
quantitatively how largely the non-Gaussianity from the magnetic field is allowed by ob- 
servations, we decide to adopt the conservative bound that the curvature perturbation 
from the magnetic field is smaller than 10~ 5 . This puts a non-trivial lower bound on e as 



e > 10 5 x 



$Pem(t 



end J 



Pinf 



As we will check later, inclusion of the second term of Eq. (fT7|) does not alter this bound. 
If one resorts to the inflationary magnetogenesis to explain the observed magnetic fields 
on Mpc scales, such inflationary model needs to satisfy this bound. 

In case of the standard single field inflation model, the generic bound (ITS]) can be 
rewritten into the bound on combination of Hubble parameter during inflation H m i and 
the so-called reheating parameter _R ra d defined by, 



R 



rad 



^end / Pend 



^reh \ Prch 



1/1 



cxp 



AN 

— 1 



-l + 3w) 



(19) 



where AN is the number of e-fold between the end of inflation and the time of reheating 
and w = P/p is the equation of state parameter during that era. This parameter was 
introduced and used in [35 |B6l - H8] and roughly measures the degree of deviation from 
the radiation dominated expansion of and duration of that era. Basic requirement that 
the inflation energy scale is less than (10 _5 Mp) 4 and the reheating occurs well before the 
Big-Bang Nucleosynthesis [H] restricts the allowed range of R ia d as —35 < logi? ra d < 12. 
Using the standard formula for the curvature perturbation [15] and observational data [39] : 



Vr 



87r 2 Mj,e 



2.4 x 10" 



(20) 



we can replace e in Eq. (118]) by H in {. Assuming no entropy is produced after reheating 
and using the reheating parameter, Eq. (ITS]) can then be written as 



M, 



nf J R r 2 ad >0.14x 



1/6 / \ 1/2 
9*,vch \ / PB,0 



9*,0 



P 7 ,0 



9.0 x 10" 



13 ( 9*,ieh 

200" 



1/6 



B 



10" 17 G 



(21) 
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where g* represents the relativistic degrees of freedom and p 7 ,o(p_B,o) is the energy density 
of radiation(magnetic field) today. This basically puts a lower bound on the energy scale 
of inflation. In particular, if the inflation end is immediately followed by reheating or 
w — 3, in which case we have i? ra d = 1, the bound becomes 

H M > 2.2 x 10 6 GeV ( ^) ^ ( -—^^ . (22) 



200 / V 10 ~ 1?G . 

In the case the energy density of the magnetic field is diluted by £ < 1 compared to the 
dominant energy density due to the entropy injection after the reheating, a factor ^~ 1 / 2 
must be multiplied on the right-hand side of Eq. (12 Therefore, if future observations 
detect tensor-to-scalar ratio, which uniquely fixes H m {, and provide an decisive answer 
that the extragalactic magnetic fields of B ~ 10~ 17 G are of inflationary origin, we get a 
lower bound on £, namely, the maximum amount of entropy injection. 

So far, we have focused on the first term on the right-hand side of Eq. ( Tl4|) . As for 
the second term, it does not dominate over the first term. To see this, let us first notice 
that the integrand of the first term contains p + P in the denominator which enhances the 
integrand by e _1 while the second term does not have such an enhancement factor. As a 
result, the contribution from inflationary era to the second term is suppressed by the slow- 
roll parameter e compared to the first term. The second term receives an enhancement 
oc log a in the radiation dominated era, but this boosts the second term by at most 
^Vrad < 100, the number of e-folds measured from the time of reheating. The second term 
can marginally reach the first term if e is as large as (9(0.01) but becomes much smaller 
for smaller value of e. Therefore, apart from 0(1) modification that may arise for inflation 
models giving large e, inclusion of the second term does not alter our bound ( fl8l) . 

Finally, before closing this section, we provide expressions for other perturbation vari- 
ables on the uniform density slice. The Einstein equations and conservation laws can be 
recast into equations that give A, a g and v in terms of £. Using the solution (fT4|) for £, 
on super-horizon scales, they are given by 

A—^, (23, 

k 2 a g = — I dt x a 3 (ti)II(ti), (24) 

ft 



a:' 



3a5(p + P) J u 
These are at most the order of ( or smaller than ( 



J dt x a 3 (ti)n(tx). (25) 



3 Calculation in the Newtonian gauge 

In this section, we again consider the evolution of metric perturbations in the Newtonian 
gauge which is widely used in the literature including j5(JJ[5T] which obtained a physically 
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incorrect result. In the Newtonian gauge, the perturbed metric is given by 



[1 + 2§)dt 2 + a 2 (t)(l - 2V)6 ij dx i dx j . 



(26) 



There is no gauge degree of freedom left in this coordinate system. In this gauge, we find 
that the evolution equation of \I/ can be written as 



3P 



P P 



4 H — - + 3H — * + ci — ¥ = S. 



P 



P P 



(27) 



where c 2 is the sound speed of the dominant component (e. g. , (? a = K for the perfect 
fluid with barotropic equation of state and c 2 = 1 for the canonical scalar field) and S is 
defined by 



,5' 



8irGa 2 H 
pk 2 



4ttG ( 5P Tcl - -n 



(28) 



Following the standard procedure (e. g. , the one given in [15]), the inhomogeneous solution 
in the superhorizon regime (k -C aH) is found to be 



a(t) 



dhS(t 



i. 



+ / dt x 



t* 



p(tx) + P(hj JU 

ajh) (pjh) + Pjh)) 
P(ti) 



h dh aih){p{t 2 ) + P{t 2 )) 



p(h 



, dtlS{tl) pMTpK) 



(29) 



Substituting Eq. (1281) and performing integration by parts, this can be further recast into 
another form which makes its time dependence more transparent: 



8nGH(t) 
a(t)k 2 
, H(t) 



a(t) 



dt\ 



dt t a 3 (ti)n(t! 

H{h) 
p(ti)+P(h) 



a 2 (Qp(QU(U] 
p{U) + P{Q 
2 



P(t 



SPrcl(t 



n(ti 



a(t) 
W) 



dt 2 a(t 2 



(30) 



This is our final form of \& induced by magnetic field. 

Now, let us first consider the behavior of this solution during inflationary era. Just 
for simplicity, we assume the single field inflation model and treat II, p, P and hence e 
too as very slowly changing variables. Then, the first term of Eq. (130]) gives the dominant 
contribution. Neglecting the other terms, Eq. (130]) during inflation becomes 



MGa 2 ^ 

* - — r^n 

3k 2 



aH\ z U 



k 



P 



(31) 



This shows that the curvature perturbation in the Newtonian gauge grows very rapidly 
in proportion to oc a 2 . Since II has the same order of magnitude as the energy density 
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of the magnetic field, \P is suppressed by a small number, which is the fraction of the 
magnetic field energy density to the total one. But it is also enhanced by another big 
number which is square of the ratio of the wavelength of the mode to the horizon radius. 
As a demonstration, if we approximate U/p at the end of inflation by the ratio of magnetic 
field energy density to the CMB energy density today and assume modes corresponding 
to 1 Mpc today have expanded 50 e- folds after the horizon exit, \1/ at the end of inflation 
becomes as large as 3 x 10 20 for B = 10~ 17 G today. This rough estimate is enough to see 
that \I/ induced by inflationary magnetic field corresponding to 10~ 17 G today becomes 
quite large at the end of inflation. 

Let us next consider the time evolution of ^> in the radiation dominated era, assuming 
that reheating completes instantly just after inflation. Using the conservation law Ila 4 = 
const, in the radiation dominated Universe, we have 



Apart from the numerical factor of 0(1), this takes the same form as that in the infla- 
tionary era (see Eq. (13TT)). However, their time dependence is quite different. Indeed, in 
the radiation dominated Universe, ^first decays in proportion to a~ 2 , which significantly 
reduces the amplitude of ^a TS t that was quite large at the time of reheating. In particular, 
at the time of horizon reentry, ^a Ist becomes as small as 11/ p. For the magnetic field 
corresponding to B = 10~ 17 G today, this becomes U/p ~ 10 -23 , which is negligibly small. 
Therefore, contributions from the remaining terms may eventually dominate over \1/ at 
late time, especially at the time of horizon reentry. Given that the anisotropic stress and 
the energy density of magnetic field are of the same order of magnitude, we can show that 
magnitude of \1/ at the time of horizon reentry is roughly given by Eq. f fT6|) . Since what we 
observe is \1/ after the horizon reentry, the observational bound \1/ < 10 -5 must be imposed 
at the time of horizon reentry, which yields the same bound as Eq. ( TT8"j) . Therefore, as it 
should be, we have confirmed that consideration in the Newtonian gauge yields the same 
bound as the one in the uniform density slice. 

The feature \1/ oc Ila 2 //;; 2 during inflation was also found in [501 EI] - It was then 
claimed in [50|[5Tj that the large amplitude of \l/ at the end of inflation generically enters 
a constant mode in the radiation dominated era, making the Universe eventually highly 
inhomogeneous and spoiling the success of standard cosmology based on the FLRW metric. 
This observation was derived by using matching conditions formulated in [52] under an 
assumption of instant transition from the end of inflation to reheating. However, as we 
have shown above, we do not find such contamination to the constant mode. Instead, we 
observe that the growing magnetic mode ^/ oc Ua 2 /k 2 , which was significantly enhanced 




(32) 



Therefore, the contribution from the first term to \l/ in Eq. (I3"U|) becomes 




(33) 
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by the end of inflation, is taken over only by the decaying mode in the radiation dominated 
era. Therefore, the breaking of the FLRW metric in the radiation dominated era is not a 
generic consequence of the inflationary magnetogenesis. 

We can also rederive the time evolutionary behavior of \& by connecting it with per- 
turbations evaluated in the uniform density slice using the gauge transformation, which 
enables us to understand intuitively the origin of the rapid growth \t> oc a 2 during inflation. 
Noting that the shear a g is identically zero in the Newtonian gauge, time translation T 
connecting the uniform density slice and the Newtonian gauge is uniquely determined by 
<y g evaluated in the uniform density slice (see Eq. fl24p ); 



8ttG 



T = J dh a 3 (ti)n(t!). (34) 

Assuming II is almost constant during inflation, we see that T grows in proportion to 
a 2 . This means that t = const, hypersurface in the uniform density slice and that in the 
Newtonian gauge deviate more and more with time. Using the gauge transformation rule 
for if), the curvature perturbation \l/ is then given by 

, 8irGUa 2 

* = ~ c + —F- (35) 

The second term grows rapidly in time and eventually dominates over The Universe 
becomes apparently inhomogeneous with time in the Newtonian gauge, but this is simply 
because the time translation (|34|) becomes very large. From a geometric point of view, 
the rapid growth of T reflects mismatch between spacetime having anisotropic stress even 
on large scales and the Newtonian gauge in which the spacial metric looks isotropic. 
Therefore, the rapid growth of inhomogeneity does not occur on the slicing such as uniform 
density slicing where the shear a g does not vanish. 

In the radiation dominated era, calculation of Eq. (134)) shows that T becomes constant 
in time. Since H decays in proportion to a~ 2 , the gauge transformation rule for the 
curvature perturbation tells us that difference between \l/ and ( becomes smaller and 
smaller with time. In particular, \1/ and ( become of the same order of magnitude at the 
time of horizon reentry. 



4 Summary 

We have studied time evolution of metric perturbation induced by magnetic field that is 
supposed to be produced during inflation. Treating the energy-momentum tensor of the 
magnetic field as first order perturbation, we solved the linearized Einstein equations on 
FLRW background spacetime. The resultant metric perturbation is at least of the order 
of 5ps/(epi n f), where 5pB is the typical amplitude of the magnetic field at the end of 
inflation, e is the slow-roll parameter and p in f is the total energy density during inflation. 
This perturbation must not exceed the observed amplitude ~ 10~ 5 , from which we could 
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obtain a generic bound on inflation models, namely, inequalities (I18p and (l2Tj) . Any 
inflationary model that achieves magnetogenesis must satisfy these constraints. 

We also performed perturbation analysis in the Newtonian gauge in which perturba- 
tions were found to rapidly grow during inflation on super-horizon scales. This weird 
behavior is simply an artifact of mismatch of the Newtonian gauge which uses isotropic 
coordinate and the anisotropic stress of magnetic field that persists even on superhorizon 
scales. Thus we do not have to worry about the breakdown of the FLRW background 
contrary to the claim of [50],[5T]. The apparently enhanced perturbation by the end of 
inflation in this gauge starts to attenuate in the subsequent radiation dominate Universe. 
The requirement that the perturbation must be smaller than ~ 10~ 5 at the time of horizon 
reentry leads to the same bound as above. 
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